83 


OG 7.2-18 


The Role of Cosmic Rays and Alftfen Waves 
in the Structure of the Galactic Halo. 

M.K. Dougherty and J.F. McKenzie (University 
of Natal, Durban, R.S.A.) and N.J. Westergaard 
(Danish Space Research Institute, Lyngby, Denmark) 


ABSTRACT 

The effect that cosmic rays and the Alftfen waves they generate 
have on the structure of the plasma distribution perpendicular 
to the galactic disk is examined. It is shown that the plasma 
distribution exhibits two length scales and the predicted values 
of gas density far from the galactic plane indicate that models 
involving hydrostatic equilibrium should be replaced by those 
allowing for a galactic wind. 

1. Introduction. It is now recognised [1,2,3] that since the energy 
densities of the interstellar gas, magnetic field and cosmic rays are 
of the same order (~ 1 eV cm -3 ) in the interstellar medium, the coupling 
between these components, which must be taken into account in a self- 
consistent fashion, leads to important dynamical effects. Here we 
examine the role of cosmic rays and the Alftfen waves they generate on 
the distribution of the plasma perpendicular to the galactic plane. 


2. General Equations . Hydrostatic equilibrium perpendicular (i.e. 
along z) requires d 

-pg, = di(Pg + Pm + Pw + *c> (1) 

in which p = interstellar gas density with pressure pg, p c = cosmic 
ray pressure, p m = magnetic pressure, p w = Alftfen wave pressure and 
g z is the z component of the galactic gravitational acceleration given 
by -3ct>(r,z) where <^(r,z) is an axisymmetric galactic gravitational 
9z 

potential such as [3]^ G M- 

4>(r,*) = [r 2 + (a . + + b. 2 )*) 2 ]* (2) 

a,, b,; a 2 , b 2 are parameters for the bulge and disk components of 
the galaxy, Mj is the corresponding mass parameter, G the gravitation 
constant, and r the radial and z the axial coordinate. 

Treating the cosmic rays hydrodynamically, that is, in terms of their 
pressure and effective diffusion constant k , the effects of convection 
and diffusion are described by the cosmic ray transport equations 1 4 J 


£c • V^SL + 0 


dz 


dz 


Y C 


F c = (Y c -1) 


VP, 


(Y c -D 


dz 


( 3 ) 


( 4 ) 


Y c = adiabatic index, F c = cosmic ray enthalpy flux, V = Alftfen wave 
velocity assumed a i. The wave energy exchange equation which 

describes how the wave pressure is generated by the cosmic ray pressure 


gradient [4,5] is 


d dp c 

4 < v2 Pw> ■ - v 5^ 


-L 


( 5 ) 
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In the absence of additional energy losses or gains, Q = L = 0 (3) 141 

and (5) give the total energy flux conservation law as ’ 


(V 1 ) 


vp. 


K 

(V 1 ) 


d Pc 

+ V2 P W = constant 


( 6 ) 


Results for the cases of strong and weak scattering . 

(a) Strong scattering (k + 0) if the cosmic ray diffusive flux is 
small compared with the cosmic ray enthalpy flux the cosmic ray trans- 
port equation (3) integrates immediately to yield the adiabatic form 


Yc _ 


or 


P c v 
P c a P^r 


= constant 


( = P 




Y r = Vs) 


(7a) 

(7b) 


e ^ r 9Y 1 e q uat ^n (6) then tells us how the wave pressure p is 
related to the plasma density, namely H w ,b 


p w p wo ^ Vp 0 )^ + Ppo 2 (y^i)( P /p 0 ) J (i"( P /P 0 ) 2 ] 


(8) 


where the subscript "0" refers to the value of a quantity at some 
reference level, eg. at the galactic plane. It is often assumed that the 

magnetic pressure is proportional to the 
plasma pressure [3,6] which in turn we 
shall assume is related to the density by 
some polytropic law of the form 


l-O-i 


. 75-1 



> 

• go 


“ ( P / p o) Yg 


( 9 ) 


Under these assumptions the magnetic 
pressure may be simply incorporated into 
the gas pressure in the form of an 
effective sound speed also given by 


. 25 


2 _ 


Y q( p qo + p mo) 


(10) 


"'so P 0 

The above "equations of state" for p , 
. Pw and p_ enable us to integrate the 
* hydrostatic pressure balance equatii 
5^ immediately to yield 


9 s 

on 


Figure 1 : 

Distributions for r = 5k pc, 
= 13,6 kms-T, 

3.8 kms " \ 


x _ _ 

CO 


so 


Ve = 500 kms -1 ; where z' = z /b„ 
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*(r,o) - *(>-,2) - n-( c /p 0 ) Yg1 ] ♦ a^KVp ) 4 -1 ] 

* !col( p o/ p ) i n- ( p / Po ) 

(Y..-1) p 


( 11 ) 


Thus for a given shape for the galactic gravitational potential (such as 
equation (2) ), equation (11) determines the plasma distribution as a 
function of r and z. Figure 1 shows the distributions of P/p 0 , Pc /P co 
and Pw/Pco perpendicular to the galactic disk. It is interesting to 
note the qualitative agreement between the theory and observations [6,7] 
denoted by crosses. The main point to note is that the model predicts 
in a natural fashion that the plasma distribution will exhibit two length 
scales; a fast decay near the galactic plane (determined by the plasma and 
magnetic pressure)followed by a much slower delay determined by the 
self-excited Alftfen wave pressure. 


( b) Weak Scattering (k large) If we assume that the cosmic ray 
enthalpy flux is negligible compared to the diffusive flux in equations 
(3) and (4) then the wave energy equation may be written 

d(V2 p ) - (Y^O) VF 0 - (12) 

az k 

where F 0 is the (constant) diffusive cosmic ray flux. 

This equation (with some assumptions for the behaviour of k) along with 
hydrostatic pressure balance generalizes the galactic halo model of 
Ghosh and Ptuskin (GP) to include Alftfen waves generated by the cosmic 
ray pressure gradient. The GP model and this generalization suffer from 
the defect that sufficiently far from the galactic plane the density 
(and therefore the gas and magnetic pressures) and wave pressures begin 
to increase and the cosmic ray pressure goes negative. This defect in 
the model arises because it is not consistent to simply assume that the 
diffusive flux is constant since in fact overall energy flux must be 
conserved (equation (6)). Thus the wave energy flux can only grow at the 
expense of the diffusive flux with the implication that the system will 
eventually drive itself into the strong scattering mode discussed above. We 
therefore expect that the results for k 40 will be similar to those for 
k = 0 with some modifications which will maximize around the point of 
inflection of the distribution of p c with z. 
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Conclusions 

The strong scattering limit predicts two length scales describing the 
gas distribution; the larger one which prevails above the plane arises 
because the cosmic rays and the waves they generate lift the gas up and 
stretch it out far beyond the galactic plane, asymptotically approaching 
a constant value pco ~ p 0 (y c - 1) 2 / (Ve/acoT (where V * is the galactic 
escape speed ). Although Poo/Po « 1 this asymptotic value for the 
density is nevertheless much greater than typical intergalactic densities 
(2.xlO" 8 xPo) with the implication that the hydrostatic model should be 
replaced with a galactic wind model. 
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